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Abstract. We present the detailed study of the stellar and gaseous kine- 
matics of the luminous early-type galaxy NGC 524 derived from the long-slit 
spectroscopic observations obtained with the Russian 6-m telescope and the 
IFU data from the SAURON survey. The stellar line-of-sight velocity distri- 
bution (LOSVD) of NGC 524 exhibits strong asymmetry. We performed the 
comprehensive analysis of the LOSVD using two complementary approaches 
implemented on top of the nbursts full spectral fitting technique, (a) a non- 
parametric LOSVD recovery and (b) a parametric recovery of two Gaussian 
kinematical components having different stellar populations. We discuss the 
origin of the complex stellar LOSVD of NGC 524. 
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1. INTRODUCTION 

NGC 524 is a luminous (Mb = —21.7 mag) lenticular galaxy settled in the 
centre of a rich group (Garcia 1993) containing a X-ray hot gas component with 
a slightly lopsided distribution with respect to NGC 524 (Mulchaey et al. 2003). 
The galaxy demonstrates nearly circular isophotes, e < 0.05 (e.g. Magrelli et al. 
1992), so from the photometric point of view it looks face-on. However, kinematic 
measurements revealed quite fast rotation of the galaxy (Sil'chenko 2000, Simien 
& Prugniel 2000, Emsellem et al. 2004) that is inconsistent with the photometric 
inclination of < 18 deg. Moreover, NGC 524 possesses a roundish disc of ionised 
gas extending up to 4 kpc (25 arcsec) from the centre (Macchetto et al. 1996) 
rotating even faster than the stars. By comparing the rotation of the ionised gas 
and the stars in the centre of NGC 524, Sil'chenko (2000) suggested that its gaseous 
disc was inclined with respect to the stellar one. The recent structural analysis of 
NGC 524 (Sil'chenko 2009) reveals that the bulge of this giant lenticular galaxy is 
rather modest and is seen only inside R 10 arcsec, while the general structure is 
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Fig. 1. Left: The positions of the SCORPIO slits and the SAURON mosaic field of 
view are overlapped onto the y-band image of NGC 524 obtained at CFHT. Right: The 
Gauss-Hermite coefficients characterizing the asymmetry of the stellar LOSVD. 



dominated by two exponential stellar discs with the scalelengths of 0.9 and 3 kpc 
respectively. 



2. OBSERVATIONS AND DATA REDUCTION 

We used two datasets derived from long-slit and integral-field spectroscopy. 

The long-slit spectroscopic observations were obtained with the SCORPIO 
universal spectrograph (Afanasiev & Moiseev 2005) at the prime focus of the 
Russian 6-m BTA telescope of the Special Astrophysical Observatory. We observed 
NGC 524 in six slit positions (P.A. = 137,217,240,295,307,317 deg, see Fig. 1 
left) going through the galaxy centre. But here we present only the long-slit 
spectrum obtained along the kinematical major axis (P.A. = 217 deg). We used 
the "green" (480-550 nm) and "red" (610-710 nm) spectral setups covering strong 
stellar absorption features as well as the emission lines of the ionised gas, Hf3, 
[Om] (A = 4959,5007 A), Ha, [Nn] (A = 6548,6583 A), [Sn] (A = 6716,6731 A) 
providing the spectral resolution of 0.22 nm and 0.31 nm correspondingly. 

Data reduction for the long-slit spectral data of NGC 524 was identical to that 
of the lenticular galaxy NGC 7743 presented in Katkov et al. (2011a). Briefly, 
the data reduction steps included: bias subtraction, flat fielding, cosmic ray hit 
removal, building the wavelength solution using arc-line spectra, constructing the 
spectral line spread function (LSF) variation model using twilight spectrum and 
night sky spectrum subtraction taking into account the LSF variation along and 
across the wavelength direction (Katkov & Chilingarian 2010), and adaptive bin- 
ning along the slit in order to achieve the minimal value of the signal-to-noise ratio 
S/N = 20 per spatial bin. 

We also used the data obtained with the integral-field spectrograph SAURON 
(Bacon et al. 2001) at the 4.2-m William Herschel Telescope. NGC 524 was 
observed in three positions of the SAURON lenslet array with 0.94 arcsec sampling 
(see Fig 1 left). The covered spectral range was 480-540 nm with the spectral 
resolution 0.48 nm. For our analysis we used the science-ready data cube kindly 
provided by E. Emsellem that we binned adaptively using Voronoi tessellation 
(Cappellari & Copin 2003) to the minimal S/N ratio of 100 per bin. 
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Fig. 2. Left: The position- velocity diagram. White line presents a radial velocity 
profile obtained using one-component model. Red and blue crosses - "bulge" and "disc" 
components in the two-component model. Right: The LOSVD at radii at -25,0,25 arcsec. 

3. DATA ANALYSIS 

3.1. SSP-EQUIVALENT PARAMETERS AND EMISSION LINE KINEMATICS 

We derived the parameters of internal kinematics and stellar populations of 
NGC 524 by fitting high-resolution PEGASE.HR (Le Borgne et al. 2004) simple 
stellar population (SSP) models against our spectra using the NBURSTS full spectral 
fitting technique (Chilingarian et al. 2007a, b). We determined SSP-equivalent ages 
T and metallicities [Z/H] of the stellar population as well as the stellar kinematics 
using the Gauss-Hermite parametrization up-to the 4th moment, i.e. v, a, hs and 
hi (van der Marel & Franx 1993). The derived parametric LOSVD exhibits a 
strong asymmetry leading to the non-physical values of ft, 3 and ft.4 (see Fig 1. 
right) corresponding to significantly negative LOSVD "wings" . 

The emission line spectrum at every spatial bin was obtained by the subtraction 
of the stellar contribution (i.e. the best-fitting model) from the observed spectrum. 
Then we fitted it with pure Gaussians convolved with the LSF in order to determine 
the kinematics of the ionised gas and emission line ratios. 

3.2. NON-PARAMETRIC LOSVD 

We propose a non-parametric recovery technique based on the full spectral 
fitting requiring no a priory LOSVD knowledge. The logarithmically rebinned 
model spectrum, J-(X) is the convolution of the assumed normalized LOSVD, 
C(v), with the rest-frame SSP model, T r {w): 



where w = ln(A), u = ln(l + v/c). We used the output SSP model of the NBURSTS 
fitting as a template spectrum J>. This equation can be considered as a linear 
inverse problem whose solution is very sensitive to the noise in the data. Hence, 
we chose to regularize the problem be requiring the LOSVD to be smooth. In 
order to do so, we use the cubic penalization V(C) = C T ■ T> T ■ T> ■ C, where T> - 
is the third-order difference operator. The function to be minimized is given by 
X 2 + XV(C). For discussion on the choice of A see Press et al. (2007). 

Using this technique, we confirmed a strong asymmetry of the NGC 524 LOSVD. 
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Fig. 3. SCORPIO long-slit spectroscopic data. First row - one-component spectra 
fitting with S/N ratio equal to 20 in the green spectral domain. Second row - two 
component spectra fitting in green domain (S/N = 40). Black points correspond to 
"disc" component, blue points to "bulge" one, red line shows Ha kinematics. 

Fig. 2 displays the result of the LOSVD recovery for the long-slit data along the 
kinematical major axis as a position- velocity diagram (left). The LOSVD profiles 
at radii —25, and 25 arcsec are shown to the right. 

3.3. TWO-COMPONENT PARAMETRIC LOSVD RECOVERY 

Another approach we use is a full spectral fitting using a two-component model 
where different stellar population components have two different pure Gaussian 
LOSVDs. An optimal template is represented by the linear combination of two 
SSPs each convolved with its own LOSVD, hence the x 2 value is computed as 
follows: 

2 v (Fi - Py ■ Eg % ■ S(T 3 ,Z 3 ) ® C(vj, a,)) 2 
X 2s 5Ff 

where £(v, a) - pure Gaussian LOSVD; F{ and 5Fi are observed flux and its 
uncertainty; S(Tj,Zj) is the flux from the j-th synthetic spectrum of SSP with 
given age Tj and metallicity Zj ; P p is multiplicative Legendre polynomials of order 
p for correcting the continuum which determined at each step of minimization 
loop by solving the linear least-square problem; kj is the j-th component weight 
(normally found by the linear minimization). The important point in this study 
is that we fixed the relative SSP contributions kj to the values derived from the 
light profile decomposition. 

4. RESULTS AND DISCUSSION 

4.1. KINEMATICS 

At all radii the LOSVD of NGC 524 clearly demonstrates the presence of at 
least two components clearly visible in Figs. 1-2 confirmed independently by the 
non-parametric and parametric LOSVD reconstruction on SAURON and SCOR- 
PIO datasets. The inner region of the galaxy that corresponds to the small ex- 
ponential pseudobulge identified by Sil'chenko (2009) is very hot in a dynamical 
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Fig. 4. SAURON integral-field spectroscopic data. From left to right: velocity, veloc- 
ity dispersion, SSP-equvalent ages and metallicities. First row correspond to one com- 
ponent fitting spectra. Second and third rows present a "bulge" and "disc" parameter 
maps for two component spectra decomposition. 

sense with the velocity dispersion exceeding 300 km s _1 . 

This component disappears at about R « 20 arcsec, where we see a sharp 
drop in the velocity dispersion and radial velocities increasing outwards suggest- 
ing the counter-rotation with respect to the main disc (see blue profiles in the 
bottom panels of Fig. 3). At large radii (R > 20 arcsec), the LOSVD recovered 
non-parametrically also becomes strongly bi-modal with the secondary component 
clearly visible in Fig. 1. 

4.2. STELLAR POPULATIONS AND EMISSION LINE RATIOS 

The radial variations of the stellar population parameters derived from the SSP 
fitting of the long-slit data are shown in Fig. 3 (right-top). The mean stellar age 
is very old (T > 15 Gyr) being at the limit of the SSP model grid in the central 
20 arcsec. The metallicity in the centre is close to the Solar value decreasing down 
to [Z/H] sa —0.4 dex in the outer disc. The results of spectrum fitting using a two- 
component model are presented in Fig. 3 (right-bottom). "Bulge" (blue points) 
and "disc" (black points) components do not have significant differences of the age 
radial profile being both very old with the "bulge" reaching the limiting age in 
the models. The "disc" metallicity decreases from [Z/H] ss 0.1 dex in the centre 
to -0.2 dex at large radii. Our analysis of the virtually noiseless (S/N > 100) 
SAURON dataset (see Fig. 4, right columns) agree very well with the long-slit 
data. The mean stellar parameters of the "bulge" component are T > 15 Gyr 
and [Z/H] « —0.1 dex, while for the "disc" component they are T w 14 Gyr and 
[Z/H] « 0.1 dex. 
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The analysis of emission line ratios using the classical diagnostics [Om]/H/3 
vs [Nil] /Ha (Baldwin et al. 1981) rules out the ongoing star formation in the 
disc of NGC 524 leaving the space for the Shockwave excitation or the nuclear 
activity. However, the latter mechanism cannot explain the large spatial extent of 
the emission line region. 

4.3. DISCUSSION 

A very good agreement between the dynamically cold "disc" component of the 
stellar LOSVD and the ionised gas kinematics suggest that the gas is rotating in 
the plane of the main stellar disc. A small discrepancy of the rotation (20 km s _1 ) 
can be explained by the asymmetric drift and well corresponds to the observed 
stellar velocity dispersion of the disc component of ~ 100 km s _1 . Old stellar 
population and the emission-line diagnostics rule out both, recent and ongoing 
star formation in the gaseous disc. The second kinematical component (blue pro- 
files in Fig. 3) probably corresponds to two different structures at different radii. 
The dynamically hot inner part without much rotation is a manifestation of the 
compact central pseudo-bulge, while at R > 20 arcsec we see the presence of a 
counter-rotating disc component that is supported by the drops in the velocity 
dispersion (300 to <100 km s _1 ) and metallicity (0.0 to —0.3 dex) profiles. 

The origin of NGC 524 has to be investigated in detail using state-of-art nu- 
merical simulations. Right now we can speculate about its evolution based on the 
observational results we have. NGC 524 might have originated from the face-on 
collision of two initially counter-rotating co-planar giant disc galaxies. The metal- 
licity difference between the two components suggests the merger mass ratio of 
about 1:4. The gas in the main stellar disc of NGC 524 might have survived from 
the original galaxy or collected later from mergers with low-mass satellites or from 
the accretion from the cosmic filaments. However, its surface density is still below 
the threshold and therefore it prevents the start of the star formation. 
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